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Abstract 

Quantum analogues of the transient and steady-state fluctuation theorems are in- 
vestigated for a harmonic oscillator linearly coupled with a harmonic reservoir. The 
probability distribution for the work done externally is derived and the following 
facts are shown: (i) In the transient fluctuation theorem, there appears a quantum 
correction of order h 2 . (ii) In the steady-state fluctuation theorem, the existence of 
a quantum correction depends on the way of driving. In the uniformly dragged case, 
the classical formula holds, while, in the periodically driven case, there appears a 
correction of order h 2 . 
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1 INTRODUCTION 

In the linear nonequilibrium regime, statistical properties of fluctuations are 
known to characterize nonequlibrium states via exact identities such as the 
Green-Kubo formula and Fluctuation-Dissipation Theorems. On the contrary, 
general properties of the far-from-equilibrium fluctuations have not been well 
understood. Recently, two new identities valid even far from equilibrium were 
found: the fluctuation theorem (FT) first found by Evans, Cohen and Morriss[l] 
and the nonequilibrium free-energy equality given by Jarzynski[2]. 

The fluctuation theorem addresses the symmetry of entropy production or 
work externally done and several different versions are known. Nevertheless, 
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it takes the following general form: 

(limUlog -gg^U ,, (1) 
\ T ^+ooJ T 6 Prob(-cr) ' V ; 

where Prob(a) stands for the probability of observing an average entropy 
production rate (or an average power applied externally), a, during the time 
interval r. The r — > +00 limit may or may not be taken depending on the 
situations. 

FT with t — > +00 limit was first found numerically by Evans, Cohen and 
Morriss[l] for nonequilibrium steady states of thermostatted systems. Then, 
it was rigorously proved by Gallavotti and Cohen [3]. This version is referred 
to as the steady-state fluctuation theorem (SSFT). On the other hand, before 
the proof by Gallavotti and Cohen, Evans and Searles[4] derived the relation 
(1) without r — > +00 limit for transient trajectories starting from initial states 
obeying the microcanonical distribution. This version is refered to as the tran- 
sient fluctuation theorem (TFT) and has been extensively studied by Evans 
et al.[5]. Moreover, it was extended to stochastically driven systems [6, 7, 8] and 
to open conservative systems [9, 10, 11]. Its relation with Jarzynski equality was 
studied as well[12] and the related topics have been extensively investigated 
(see e.g., references in Refs. [11,13]). 

Recently, a beautiful experiment by Wang et al.[14] confirmed TFT for a col- 
loidal particle kept in a uniformly moving optical trap. The dynamics of the 
colloidal particle is governed by the Langevin equation and TFT holds for 
the work done externally. FT for the Langevin equation was studied by Ma- 
zonka and Jarzynski[15] to illustrate the difference between SSFT and TFT, 
and the experiment stimulates the reinvestigations of FT for the Langevin 
equations [16, 17, 18, 19]. 

Contrary to FT for classical systems, FT for quantum systems has not been 
well understood. In order to generalize FT to quantum systems, it is necessary 
to identify entropy change or work done externally. And two procedures are 
available: 

(a) Measure energy, particle numbers etc. twice and evaluate the entropy change 
or work done as the difference of the two observed values. 

(b) Measure flows of energy, particle numbers etc. and evaluate the entropy 
change or work done as accumulated values of the flows. 

Needless to say, in classical systems, the two procedures give the same values 
of the entropy change or work done, because of the conservation of energy and 
particle numbers. However, in quantum systems, as a result of the noncom- 
mutativity of dynamical variables, the two procedures are not equivalent. 
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As far as the authors know, the procedure (a) was first suggested by Kurchan[20] 
and he showed that the fluctuation theorem (1) holds for the probability dis- 
tribution function of the so-obtained entropy change or work done. Similar 
results were discussed by Hal Tasaki[21], Callens et al.[22], and the C* gener- 
alization was given by Tasaki and Matsui[23]. 

On the other hand, the procedure (b) has not been well studied. In this pa- 
per, we investigate a quantum analogue of the fluctuation theorem, following 
the procedure (b), with respect to the work done externally for a harmonic 
oscillator linearly coupled with a harmonic reservoir. The oscillator is exter- 
nally driven in such a way that the center of the harmonic potential follows a 
given trajectory. Note that the system is a quantum analogue of the Langevin 
equation[15,16,17,18,19] which describes the experiment by Wang et al.[14]. 
Then, we have found that (i) there appears a quantum correction of order h 2 
in TFT, and that (ii) the existence of a quantum correction for SSFT depends 
on the way of driving (no correction in the uniformly dragged case and h 2 - 
correction in the periodically driven case). Moreover, the quantum correction 
for SSFT is universal in the sense that it depends only on the temperature 
and the frequency of the driving, but not on the specific features of the system 
and environment. 

The rest of this paper is arranged as follows. The model is described in Sec. 2. 
The equation of motion is solved and quantum analogues to the fluctuation 
theorems are derived in Sec. 3. The subsequent sections are devoted to the 
discussions of TFT and SSFT separately. In the last section, discussions are 
given. 



2 MODEL 



In order to discuss quantum corrections to the fluctuation theorem, we study 
an exactly solvable model of a one-dimensional harmonic oscillator linearly 
coupled with a harmonic bath[24,25]. Because a straightforward derivation of 
a quantum Langevin equation is known[25], the Gardiner's version is used. 
Then, we consider the following thought experiment: 

(i) The system is prepared to be in equilibrium with inverse temperature (3 1 . 

(ii) At time t = 0, the harmonic potential starts to move, where its center 
follows the trajectory x = f(t). 

When t < 0, the potential center is assumed to be fixed at x = and to move 
continuously (namely, /(0) = 0). The Hamiltonian which describes the time 

1 Throughout this article, we use the unit where the Boltzmann constant is unity. 
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evolution is given as 



H(t) = H + ±f(t)(-2q + f(t)) , (2) 

H o = ^ + ^ 2 + \--J rfA ((PA - n x qf + ulql) : , (3) 

where q, p, and m stand for the coordinate, momentum and mass of the 
harmonic oscillator, respectively, k is the strength of the harmonic potential, 
qx and p\ are the coordinate and momentum of a bath degree of freedom with 
frequency U\, and : • • • : is the normal product with respect to the normal 
modes: 

a x = -/±=(u, xqx + i Px ). (4) 



Here we assume that u\ runs from to +oo and that the dispersion equation 
r)(z) = mz 2 — k — d\n\ — / d\ 2 A A 2 (5) 

J J Z UJ\ 



has no real zeros, where 

V\ = v^A^A • (6) 

This condition corresponds to the case of the damped harmonic oscillator. 

As in the previous works[14,16,17,15], we investigate the statistical property 
of the work externally done during the time interval r divided by temperature 
P^ 1 : E T . If the work is started to be measured at time t, it is given by 

T 

± T = pJf(s + t)(-k(q(s + t)- f(s + t)))ds (7) 



where / stands for the time derivative of / and q(s + t) is the coordinate 
operator in the Heisenberg picture at time s + t. Then, the probability density 
7Tt(S T = A) of the work is given by 

* t {Z T = A) = (e-< la °6(& r -A)) (8) 



where (•••)= tr(e P H ° ■ ■ -)/Z stands for the thermal average and Z = tre * H ° 
is the partition function. Note that TFT refers to the symmetry of 7To(S r = A) 
while SSFT refers to that of tt +00 (S t = A). 
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3 DISTRIBUTION FUNCTION OF THE WORK 



In this section, we derive an explicit expression of the distribution function 7i t 



of the work E T . 



3.1 Evaluation o/E 7 



First we note that the unperturbed Hamiltonian H is diagonalized by the 
normal modes: 

v x ip-imuxq , f dX' ( v y a x > v x >a y \\ 
a x = a x + — 7 — ^ \ j= — + / ~7T — 7T : ( ( 9 ) 



r)+{u\) { \/2h J 2 \u x -uj X i + iQ ui\ + uiy 

where the normal mode a\ and the interaction strength v\ are given, respec- 
tively, by (4) and (6), 77+ (u>\) = r]{uj\ + iO) and the function r)(z) is defined in 
(5). Indeed, one has [a\, a x ,] = 5(X — A') and 

Hq = J dXhu\a x a\ + (c — number) . (10) 
The old variables can be represented by these normal modes. For example, 



and, thus, 

Hit) - H - ^ km j d x{^ - ^L) + (c - number) (12) 



where ri^(u\) = rj{tjj\ — iO). 

From (12), one finds that the normal mode a\(t) in the Heisenberg picture 
obeys 



a x (t) = -iuj x a{t) + 13 
V2h 



which admits the solution: 

kv x 



a x (t) = a x e~^ - i ^= r ^~ ( : f(t) - f e^*-*' f(s)ds (14) 
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And we obtain the coordinate q(t) in the Heisenberg picture as 



a x 7 — ra 



+ k f dX , V ( X ,, 2 / sin(cu A ( S - (15) 

Then, by a tedious but straightforward calculation, one obtains the work op- 
erator XL 



t+T S 



+ / dX f k Y x [ ds [ ds'cos(u x (s-s'))f(s)f(s') (16) 
J uj x I r?+(cj A ) \ 2 { J 

where f(u\;t,r) is a 'partial Fourier transformation' of fit): 

t+T 

f(u x ;t,T)= I f(s)e-^ s ds. (17) 



t 



3.2 Calculation of the work distribution n t 



The distribution function 7r t (X T = A) of the work S r is obtained with the aid 
of the characteristic function 

$ t (0 = JdA 7r t (S r = A) S A = (exp (i£± T )) . (18) 



With the aid of the formula 

( e /dA(&a A +T, A a+)) = JdA^coth^) ^ 



(its proof is given in Appendix), one obtains the following expression: 

$ t (f) = e ^.--#<r ? (20) 
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where 

mtir = „*» 'Yds [ U I dX . f, l .r^ (S ~ S ' )) fW(s') (21) 



and 

= ^rJ ds J ds J dX \l(L)p - c ° th (V )/(8)/(8) - (22) 



Eq.(20) shows that the distribution 7r t is Gaussian with the mean value m tjT 
and the standard deviation <r tiT : 

^--^Tlk-t- 1 ^)- (23) 



By setting £ = 0, (23) leads to the quantum counterpart of TFT [4,5] 

1o Oe St= 1 =i?^- < 24 > 



In addition, provided that the limits m +OO T = lim^ +00 m t;T and <7 +0O)T = 
lim t ^ +00 (T t r exist, one has 



1 ( A ~ 



m_ 



\2' 



7T +00 (S r = A) = — = exp —2 ; — . (25) 



where 7r +00 (S r = A) = lim t _,. +0O 7r t (E T = A), and the quantum counterpart of 
SSFT is derived 

km - log — = hm a . 26 

t-,+co T 7T +00 (S r /r = -a) r^+oo a 



If the coefficient of A in (24) and that of a in (26) were unity, the fluctuation 
theorems would hold. However, the expressions of the mean value m ttT and 
the standard deviation a t>T seem to indicate the deviations. In the following, 
we investigate this question for TFT and SSFT separately. In the latter case, 
since the existence of the t — > +oo limits strongly depends on the behavior of 
fit), we focus on two concrete examples: the uniformly dragged case fit) = vot 
and the periodically driven case f(t) = v /QsmQt. 
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4 QUANTUM CORRECTION FOR TFT 



In this section, quantum corrections for TFT (transient fluctuation theorem) 
will be investigated. 

First we note that a\ T can be rewritten as 

el = tt¥ Ids Ids' /dA "> C ° s( ^-/' )) coth(^)/ M /( S -).(27) 



As the integrand is different from that of 2m 0jT , one has trj r ^ 2mo )T and, 
thus, TFT does not hold. 

Using the expansion coth(^| ;2L ) = + \l3hu)\ + 0(/i 3 ), the quantum cor- 
rection can be easily calculated. Indeed, one has 



*0,r 2 = ttf? ids [ds> I dX , f^^^'^ hsms') 

^pWjds j ds'J dX J ( A J* cosMa - s'))f( S )f( S ') 







+ 0(£ 4 ) (28) 
and, thus, 

log *fy T =\ = (l + ^ + 0(n 4 )) A . (29) 



where 



62 = "24^7 / dS I dS ' I dX \ V +(Z\ 2 COS( " A(S ~ S ' ))/(S)/(S,) • (30) 



o o 



To illustrate the behavior of e 2 , we consider the uniformly dragged case f(t) = 
v t where the damping is very weak: /i = J K^dX/k <^ 1. Up to the first order 
in /x, the function l/r)±(u>) can be approximated by 



(31) 



r] ± (to) m (uj - Q Q ± ij) (uj + VLq ± ij) 
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= # + / ^ I + 1 (32) 

V m J 4m [ _ Ux ^ k / m + CUx j 

1 = inl dXK l 5 ("J k / m ~ ^0 ( 33 ) 

where the symbol P in the integrand stands for Cauchy's principal part. Then, 
by a straightforward calculation, one finds that the leading term of the cor- 
rection e 2 with respect to ji is independent of time r and 

£ --^ (1 + 0( ")»' < 34 » 



5 QUANTUM CORRECTION FOR SSFT 



In this section, quantum corrections for SSFT (steady-state fluctuation the- 
orem) will be investigated. In this case, the long-time limit t — > +oo should 
be taken. As the existence and the limiting value strongly depends on the 
potential motion fit), we examine quantum corrections for two examples: the 
uniformly dragged case f{t) = v t and the periodically driven case f(t) = 
Vq/Q sin fit. 



5.1 Uniformly dragged case 

In this case f(t) = vot and the mean work m tiT and its variance of T read as 
m tT = 2(5k\l fdX , V ) ^(^(t + r/2))sin^/2 ^ 



and 



't,T 



2hk 2 (3 2 v 2 / dX 



-coth( 



a v sin- 



u x r/2 



(36) 



In the limit of t — > +oo, one finds 



, f r 2-k [3k 2 vfal sin lu x t/ 2 
m +00>T = lim m t:T = / dX- — 5{u X ) = M +00 r (37) 
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where v\ = y/tJ~\K\ is used and 

M +0 o = / dX^SM (38) 

Note that a\ OOT = a 2 T as the variance does not depend on t. 

These equations together with {sin 2 ut/2}/(u 2 t) — > nS(u)/2 (r — > +00) give 



a 2 -Khk 2 l3 2 vl r k\ (3hw x 

hl fL = oT7 / rf V M2 ^ coth(^— )5(a; A ) 



r-+oo2m +00iT 2M +oc y I ?7 + (^ A ) | 2 " 2 
nhk 2 (3 2 v 2 2 



2M +00 /3nj 77+(0) 
Therefore, SSFT holds exactly: 



dXK{5(uj x ) = 1 . (39) 



i- 1 ^ vr +00 (S r /r = a) 

hm -log — — r — a. 40 



The reason can be understood easily. Indeed, the motion of the potential is 
uniform in this case and only the zero frequency environmental modes do 
contribute to the steady-state work distribution. On the other hand, quantum 
effects appear only for the non-zero frequency modes and, thus, the quantum 
correction disappears in SSFT. In the next section, we consider a case where 
non-zero environmental modes do contribute to SSFT. 



5.2 Periodically driven case 

In this section, we consider the periodically driven case fit) = vocos(Qt). 
First we remark that, in this case, there exists no exact steady state since the 
system is periodically driven. We consider, instead, a limiting oscillatory state 
obtained by a 'stroboscopic limit', where the time t is set to t = 27rm/Q + 4>/Q 
and the limit m — > +00 is taken. The variable <p represents the phase of the 
limiting oscillation. Hereafter, we abbreviate 

^F(t)|^ f+ ,=strob. t lim F{t) (41) 
The variance of the work for t = 2nm/il + <p/Q, now reads as 



aL= ne^vi j ds j d , j dx v 2 cos(u x (s - s >)) 
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xcoth( ^ A ) cos(f2s + 0) cos(f2s' + 



(42) 



Since it does not depend on m, one has c+oo,t = strob. lim t _ >+00 a 2 T 
Then, by a tedious but straightforward calculation, one obtains 



>t,r- 



lim cr +0OT 



7rM; 2 /3 2 t; 2 coth^jp 
4 



(43) 



As in the same way, one has 



m +oc T = strob. lim m t T 

t— »+oo ' 



v 2 x cos(u x (s - s')) 



j3k 2 Vn ds ds / d\— r——j — r- t^t~ cos(f2s + 0) cos(f2s' 

J J J u x \ V+i^x) \ 



Tl(5k 2 V 2 



-oo 
2„,2 



4 n\rj+(n)\ 



cos 20 - cos(2fir + 20) 



fi 2 -w 2 



sin(2fir + 20) - sin 20 
2SL 



(44) 



and 



vr/3/c 2 ^ 
hm m +00jT /r = — - — — — 
r-»+oo 4ll\ri + (ll)\ 2 



J d\v 2 x 5(uj x - Q) . 



(45) 



Eqs.(43) and (45) give 
2m_|_ 00iT 



lim 



2 , 

tanh 



+ 00,T 



(46) 



and, thus, SSFT does not hold 

,. 1, 7r +00 (S r /r = a) 
lim -log —- r 

T^+oc r 7r +00 (S r /r = -a) 



2 /Ml 

tann— — a ^ a 



phtt 



(47) 



where the distribution function refers to the limiting oscillatory state. Note 
that the quantum correction is universal in the sense that it depend only on 
the temperature and the frequency of the driving, but not on the specific 
features of the system and environment. 



11 



6 DISCUSSIONS 



In summary, for a harmonic oscillator linearly coupled with a harmonic reser- 
voir, we have investigated quantum analogues of the transient and steady-state 
fluctuation theorems with respect to the work done externally. And we have 
shown the followings: 

(i) In the transient fluctuation theorem, there appears a quantum correction of 
order h 2 . 

(ii) In the steady-state fluctuation theorem, the existence of a quantum cor- 
rection depends on the way of driving. In the uniformly dragged case, the 
classical formula holds, while, in the periodically driven case, there appears 
a correction of order h 2 . 

We remind that the steady-state fluctuation theorem for the periodically 
driven case refers to the limiting oscillatory state, which is obtained as the 
stroboscopic limit defined in (41). 

Remarkably, the situation is in contrast to that for Kurchan's quantum fluc- 
tuation theorem, where energy is measured twice, and the entropy change and 
work done are evaluated from the difference between the two measured values. 
In this discussed previously [20, 21, 22,23], the (transient) fluctuation 

theorem holds as in the classical systems. 

Finally, we note that the quantum correction for the steady-state fluctuation 
theorem is universal in the sense that it depends only on the temperature and 
the frequency of the driving, but not on the specific features of the system 
and environment. We expect that this is also the case for general systems. 
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A Derivation of ( e I dX{ ^ +r "< ) ) = e pA«Ap C oth(^P)_ 

The desired result immediately follows from 

< e v + +fr) = e ^ coth ^ , (A.l) 

where [a, a + ] = 1 and (•••) = tr(- • • e - phja+a )/Z with Z = tr(e^^ a+a ). This 
relation can be shown as follows. 

The Baker-Hausdorff theorem gives 

With the aid of the coherent states \a) = J2™=o 7=7 l n )> which satisfy 

a I a) — a\a) , (a\a + = a* (a\ (A. 3) 

f — I a ) ( a I e -M 2 = 1 , (a\P) = e a *" , (A.4) 

J IT 

one finds 



£2 

z 



Ve-^(n I I— I a) (a | e" a V a f ^- \ f3')((3' \ n)e 

n=0 J 71 J n 

Z n J 71 J 7T n\ 



n=0 



_ ej^ r dcx r d(3' ^_ M 2_ ll3 , l 2 + £ l3 , +r)a * +aj3 ,* e -i3hu +a * l3l 

Z J IT J IT 

= e lT— [ d P' c -(l-e-P h "W\ 2 +e-l 3h "l3"- V +P't 
Z J IT 

= e ^ coth ^) , (A.5) 
which leads to the desired relation (A.l): 

< e" a+ +« a >= e^ coth ^ (A.6) 
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zplane 




Fig. B.l. Integration contour for the integral (B.3) 



B Integral containing r] + (u\) 



To proceed the calculations in Sec. 3, one has to evaluate integrals involving 
ri + {u\). In this appendix, we illustrate the way by evaluating 

J u x | r] + {uj x ) \ 2 

The key idea is to rewrite it as a contour integral and to use the relations: 



r)-(x) — r) + (x) = — ni J d\8{x — u\)v\ , (B.l) 
r 1± {-x)=7 lT {x) . (B.2) 



Then, one has 



2 

/ d\ — ^ — — = [ d\ [ dx5(x - lu x ) 

J LUx 77+ OA) \ J J 



u x \v+M\ 2 J J r]^(x)r] + (x) x 

i,i i.i, 7 i ii, 

-ax = / -——dz 



o 



-ixi r] + (x) f]-{x) x J — irir] + (x)x 



f 1 dz , 2 f ,„ 1 rie* , ? 1 Rie ie d6 

= / — — hm / dO T7T-, — — — urn 

J —irizriiz) r^oj 



zt)(z) r ^ ^ —iri re te i](re te ) R->oo7 —ni Re t9 i](Re' 19 ) 
- I • (B.3) 
where the integration contour C is shown in Fig. B.l. 
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